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Abstract 

The resonant formation of nonlinear coherent modes in trapped Bose-Einstein con- 
densates is studied. These modes represent nonground-state Bose condensates. The 
methods of describing the spectrum of the nonlinear modes are discussed. The latter 
can be created by modulating the trapping potential with a frequency tuned close to 
the transition frequency between the two chosen modes. The requirement that the 
transition amplitudes be smaller than the transition frequency implies a constraint on 
the number of particles that can be transferred to an excited mode. The resonant Bose 
condensate serves as a collective analog of a single resonant atom. Such a conden- 
sate, displaying the coherent resonance, possesses several interesting features, among 
which are: mode locking, critical dynamics, interference patterns, interference current, 
atomic squeezing, and multiparticle entanglement. 
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1 Introduction 



Dilute atomic gases, trapped and cooled down to temperatures when almost all atoms are in 
the Bose-condensed state, are described by the Gross- Pit aevskii equation (see reviews [1-3]). 
The mathematical structure of the latter is that of the nonlinear Schrodinger equation which, 
due to the presence of the trapping potential, possesses a discrete spectrum. The equilibrium 
Bose-Einstein condensate corresponds to the ground state associated with the lowest energy 
level of the spectrum. If not the ground but an excited state would be macroscopically 
populated, this would correspond to a nonground- state Bose condensate. The possibility of 
creating such a nonequilibrium condensate was advanced in Ref. [4]. This can be done by 
applying an alternating field with a frequency tuned to the transition frequency between the 
ground state and a chosen excited state. The latter states are called the nonlinear coherent 
modes and they are described by the stationary solutions to the Gross-Pitaevskii equation. 
The properties of these modes have been considered theoretically in several publications 
[4-13] and a nonlinear dipole mode was observed experimentally [14]. A known example of 
such a mode is a vortex that can be formed by means of a rotating laser spoon [15,16]. 

The aim of the present publication is threefold. First, we give a survey of the problems in 
characterizing the nonlinear coherent modes, especially in the accurate calculations of their 
spectrum. Second, we stress the analogy of the resonant Bose condensate with a resonant 
atom. And, third, we show that, because of its coherent collective nature, the resonant Bose 
condensate possesses a number of properties distinguishing it from a single resonant atom. 
We describe several interesting novel effects that can be observed in a nonequilibrium Bose 
condensate. 

2 Nonlinear Coherent Modes 

First of all, it is necessary to show how one could get an accurate description of nonlinear 
coherent modes. We assume that the system of N Bose atoms is confined in a trap, with a 
trapping potential U(r) tending to infinity at large r = |r|, 

U(r) — > oo (r — > oo) . 

Due to the presence of this confining potential, the spectrum of the stationary Gross- 
Pitaevskii equation is discrete, being defined by the eigenvalue problem 

H[<p n ] <p n (r) = E n v9 n (r) , (1) 

with the nonlinear Hamiltonian 

H[p] = - ^ + U(r) + N j $(r - r') |y,(r')| 2 dv' , (2) 

where m is atomic mass and <3>(r) is a potential of interatomic interactions. The eigenfunc- 
tion (p n (r), labelled by a multi-index n, is a nonlinear coherent mode. It is assumed to be 
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normalized to unity, ||y? n || = 1- Atomic interactions for dilute trapped gases are described 
by the Fermi contact potential 

$(r) = A s 5{v) , A s = 4nh 2 ^ , (3) 

with a s being the s-wave scattering length. The trapping potential is commonly presented 
by the harmonic oscillator form 

U^^^lrl + ^rl + ^rt) . (4) 

Here, we shall consider the single- well potential (4), though the double- well and more compli- 
cated potentials [6,11,13] are possible, including periodic potentials corresponding to optical 
lattices [17-31]. 

In the case of cylindric symmetry, when it is convenient to introduce the 

notation 

v = — , l r = J (5) 

for the frequency ratio v and oscillator length l r , which allows the usage of the dimensionless 
spatial variables 



\l r l + r y 

• _L = 

And the quantity 



g = An^N (6) 

is a dimensionless coupling parameter. Also, let us define the dimensionless Hamiltonian 
and the wave function, respectively, 

Then Eq. (2) reduces to 

^ = -^V 2 + i(ri + ,V)+,H 2 , (7) 

where 

„ 2 d 2 Id Id 2 d 2 
V = 1 1 1 . 

dr\ r± dr± r\ dip 2 dz 2 

The eigenproblem (1) takes the form Hip nm j = E nm jifj nm j, in which n — 0, 1, 2, ... is the radial 
quantum number, m — 0, ±1, ±2, ... is the azimuthal quantum number, and j = 0,1,2,... 
is the axial quantum number. 

The equilibrium Bose-Einstein condensate corresponds to the ground-state solution to 
the stationary Gross-Pitaevskii equation, when n = m = j = 0. This equation is nonlinear 
and does not allow for an exact solution. If the coupling parameter (6) were small, g — > 0, 
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perturbation theory would be admissible. And in the opposite limit of asymptotically strong 
coupling, g — > oo, the Thomas- Fermi approximation could be used (see [1-3]). Here we 
demonstrate how one can construct an accurate solution for the ground state, which is valid 
for arbitrary coupling parameters g G [0, oo). For simplicity, we consider the spherical trap, 
when v — 1. Then the ground state ipo(r) = V ; ooo( r ) depends on the spherical variable 
r = \Jr\ + z 2 . 

To find an accurate approximate solution of a nonlinear differential equation, we employ 
the self-similar approximation theory [32-36] , in the variant [36-38] designed for constructing 
self-similar crossover approximants satisfying the corresponding asymptotic conditions. To 
this end, we notice that at short radius r — > 0, when the nonlinear term in the Hamiltonian 
dominates, one has 

^ (r) ~ c + c 2 r 2 + c 4 r A (r -> 0) . (8) 

And at large r — > oo, where the harmonic term in Eq. (7) becomes dominant, the wave 
function tends to the Gaussian form 

tpo(r) ~ C exp ^— i r 2 ^ (r — > oo) . (9) 

The self-similar crossover approximant, interpolating between the asymptotic limits (8) and 
(9) is 

■0o( r ) = C exp ^— - r 2 ^ exp jar 2 exp(— 6r 2 )| . (10) 

The coefficients a and 6 are obtained by expanding function (10) in powers of r and substi- 
tuting this expansion in the stationary Gross-Pitaevskii equation, which yields 



a = l + \(gC 2 -E ) , 6 = -L (1 - 2a) (E - 1 + 2a) 
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10a v / v « - / 20a 

The ground-state energy Eq and the normalization constant C are defined by the consistency 
conditions 

4tt / ^o(r) ^o(r) r 2 air = E , An ^l(r)r 2 air = 1 . (11) 
jo Jo 

The accuracy of an approximate solution can be characterized by calculating the local residual 

R(r) = V^r(H-E )Mr), (12) 

showing the deviation of the radial wave function from the exact solution at each point r. 
The behaviour of this residual for the self-similar approximant (10) is presented in Fig. 1 
for different values of the coupling parameter (6), where it is compared with the residuals 
for the optimized Gaussian approximation and Thomas- Fermi approximation [3]. As is seen, 
the accuracy of the self-similar approximant (10) is essentially better. Figure 2 shows the 
atomic density 

n(r) = Aix^lij-) (13) 

for the self-similar approximant (10), as compared with that for the optimized Gaussian and 
Thomas-Fermi approximations. 
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The way of constructing an accurate approximation for the wave function is suitable 
for describing the ground state in a spherical trap, but it becomes rather cumbersome for 
nonspherical traps and for nonground states. Moreover, the exact shape of wave functions 
is of less importance than the related energy values. Hence it is necessary to concentrate 
attention on the calculation of the spectrum of higher nonlinear modes. For this purpose, 
the most suitable is the optimized perturbation theory [39-41]. This approach has been 
successfully employed for a variety of physical problems. A brief survey of the method and 
many citations can be found in Refs. [35,36]. Recently, this approach was used for calculating 
the critical temperature of Bose-Einstein condensation in a dilute Bose gas [42-47]. The 
main idea of the optimized perturbation theory is the introduction of control functions that 
optimize the convergence of a calculational procedure [39-41]. Control functions can be 
introduced in different ways, and different calculational algorithms may be employed. Any 
kind of perturbation theory or iterative procedure can be used. The Rayleigh-Schrodinger 
and Dalgar no- Lewis perturbation theories can be invoked for the Schrodinger equation as 
well as for Green functions [48]. Semiclassical expansions are possible [49]. One may also 
prefer to work not with the Schrodinger equation itself but with an inverse Schrodinger 
equation [50]. 

We apply the optimized perturbation theory, with the Rayleigh-Schrodinger algorithm, 
to the Gross-Pitaevskii equation (1). For the spectrum, in the first approximation, we find 

e(g) = E(g,u(g),v(g)), (14) 

with the expression 

E(g, u, v) = | (u + ^ + | + ^ + u y/v gl nmj , (15) 

in which 

p = 2n + \m\ + 1 , q = 2j + 1 
are the combinations of quantum numbers, and the integral 

j /*0O /"27T /"+00 

I nmj = — -= \ r ± dr ± dip dz \ip nmj (r ± ,(p, z)\ 4 

uyv Jo Jo J-oo 

contains the wave functions 

-2n!«M+i"| 1/2 

V>nmj{r±,(p,Z) = 



n + \m\)\ 

irmp / v \l/4 I 



U 



exp [ - - r± ) x 



\7lJ jWJ\ 

where L™(-) is a Laguerre polynomial and Hj(-) is a Hermit polynomial. The control func- 
tions u = u(g) and v = v(g) are obtained from the optimization condition 

( d d\ 
\du +Sv dv) E ^ 9, M ' ^ = ° ' 
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which results in the equations 



( .,,-l) + ^/ = o, 9 (i-4) + _!G =o, (16) 

u z J pis w q \ v z I pis Jvq 



in which 

G = 2p y/q I nmj gv . (17) 

Equations (15) to (17) define the spectrum (14) for all quantum numbers n, m, and j and 
for arbitrary values of the coupling parameter g. 

Instead of solving numerically Eqs. (14) to (17), we may derive explicit analytical ex- 
pressions for the spectrum of nonlinear modes by resorting to the technique of self-similar 
root approximants [36-38]. To this end, we find the weak-coupling expansion of spectrum 
(14) in the form 

e{g)~Y.a n G n (G -> 0) , (18) 

ra=0 

with the coefficients 

Also, we obtain the strong-coupling expansion 

e(g) ~ ]T b n (G — > oo) , (19) 

n=0 

where the coefficients are given by the equalities 

4b = 5 , 46i = 2p 2 + Q 2 , 

206 2 = -3p 4 + 2p 2 Q 2 - 2Q 4 , 296 3 = 2p 6 - p 4 Q 2 - 2p 2 Q A + 2Q 6 , 
500&4 = -44p 8 + 22p 6 Q 2 + 2p 4 Q 4 + 78p 2 Q 6 - 69Q 8 , 
125006 5 = 1122p 10 - 595/Q 2 - 70p 6 Q 4 + 440p 4 Q 6 - 3640p 2 Q 8 + 2821Q 10 , 
and the powers are 

2 2 6 

A) = - ft = - - , (32 = - - , 

(3 3 = -2 ft = - y , A = - f ■ 

Interpolating the asymptotic limits (18) and (19) by means of the self-similar root approxi- 
mants [36-38], we obtain 



el(g) = a ^...(l + A kl G) nkl + A k2 G 2 ) + . . . A kk G j . (20) 

Here, depending on the approximation order k — 1, 2, . . ., we have in the first order 

1.746928 2 
Al1 ~ 5/2 — ' nn ~ 5 ' 



in the second order 



(V+^)V6 3.051758 6 1 

A 21 = 2.533913 ^Te ' A ^ = 5 > n 2i = 7 > n 22 = 7 

a 7 a o oo 



and in the third order 

A 31 = 1.405455 125/22 2 ^ 2 , 5/66 > ^32 = 6.619620 ^ , 

5.331202 _ 6 11 2 

As3 ~ a J5/2 ' n31 ~5' n32 ~10' n33 ~15- 

The accuracy of approximants (20), for vortex nonlinear modes with n = j = and different 
to, as compared to the numerical values of Eq. (14), are illustrated in Fig. 3. As we see, the 
approximants e 2 (g) and e\{g) are quite accurate. 

The analytical expression (20) makes it possible a direct analysis of the spectrum with 
respect to varying quantum numbers and the coupling parameter. It also allows us to get a 
fast estimate of the critical parameter g c , and, hence, of the critical number of atoms above 
which the system of atoms with attractive interactions becomes unstable. For this purpose, 
we need to find the negative value of G c at which the spectrum (20) becomes complex. To 
illustrate the idea, we limit ourselves by the first approximation e\{g), which is complex for 
G < G c , with 



5/2 



G c = -0.572433 a . 
This, in view of notation (17), gives 

g ^-0M teg, 

which, because of Eq. (6), yields the estimate for the critical number of atoms 



N ^ (2 P + q u)^ 



300^^ I nmj v 



This estimate is of order of that obtained by means of numerical solution [3,4]. A more 
accurate value of iV c can be derived from the consideration of the higher-order approximants 
(20). 



3 Resonant Bose Condensate 

Suppose that a trapped atomic gas is cooled down to temperatures when almost all atoms are 
in the Bose-Einstein condensate corresponding to the ground state energy Eq. As has been 
shown [4], applying an alternating external field, with the frequency close to the transition 
frequency between the ground state and an excited mode, it is possible to transfer atoms 
from the ground to this excited mode. Thus, a nonground state condensate can be created. 
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The latter, by applying one more resonant field can be transferred to another higher mode. 
In this way, we may consider resonant transitions between two energy levels, say E\ and E 2 , 
such that Ei < E 2 , with the transition frequency 



^21 = ^ (E 2 — Ei) . 
The transition is induced by an alternating field 

V = Vi(r) cos cut + V 2 (r) sin cut , 
whose frequency is close to frequency (21), so that the resonance condition 



«1, 



uo — uo 2 i 



'21) 



(22) 



(23) 



holds. This situation is similar to the resonant interaction of electromagnetic field with 
a single atom, though a Bose condensate is a many-particle system and the transition is 
between collective levels. 

In the presence of the alternating field (22), we have to deal with the time-dependent 
Gross-Pitaevskii equation 



ih^-ip{v,t) = (H[ip\ + V)<p(T,t) , 



dt 



(24) 



with the nonlinear Hamiltonian (2). One may look for the solution to this equation in the 
form of the expansion 



<p(r,t) =^2cn(t) (f n (r) exp (- ^ E n tj 



(25) 



over the nonlinear coherent modes y? re (r) defined by the eigenvalue problem (1). Note that, 
in general, one could look for a solution expanding the latter over an arbitrary basis [51]. 
However, the set {tp n (r)} of the nonlinear coherent modes stands out as the sole natural 
basis, being physically distinguished from all other formally admissible expansion sets. 
There are two types of transition amplitudes in the problem, 



N 



\<Pr. 



2|y?n(r) 



m. 



dr , 



Pmn^lJ </4(r) [Vi(r) - iV 2 {v)\ <p n (r) dv 



(26) 



the former being due to atomic interactions (3) and the latter, to the action of the modulating 
field (22). In analogy with the resonance in an atom, we need here that the transition 
amplitudes (26) be smaller than the transition frequency (21), 



a 12 

^21 



«1, 



«21 
^21 



«1, 
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^21 



«1, 



(27) 
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in order that the transition frequency, and the resonance as such, be well defined. These 
conditions lead to the restriction 



< (V + Q>)«< 

lApy/q I nm j 



on the system parameters, which can be reformulated as the limitation 

(2p 2 + Q 2 ) 5 / 4 



N < 



hQnvp^/q I nmj 



(29) 



on the number of particles that can be resonantly transferred to a chosen mode. Note that 
the limiting value in the right-hand side of Eq. (29) is close to iV c . 

Under conditions (27), the Gross-Pitaevskii equation (24) reduces to the system of equa- 
tions 

1 ~dt = ai2 ' C2 ' Cl 2 ^ 12C2e ' 1 ~dt = tt21 ' Cl l ° 2 + 2 l2 ° ie ^ ' 

for the fractional population amplitudes. This reduction to an effective two-mode system has 
become possible due to the resonance condition (23) and inequalities (27). Such a situation 
is analogous to a resonant atom whose description also reduces to a two-level system. One 
more requirement is that the energy levels would not be equidistant [4]. This is always valid 
for electronic levels in an atom. Fortunately, owing to the nonlinear term in the Gross- 
Pitaevskii equation, this is also true for the spectrum of the nonlinear modes. There is as 
well another possibility of shifting the required levels by a slight modification of the trapping 
potential [52]. 

Equations (30), constituting a system of four equations for real functions, allow a further 
simplification. For this purpose, we may write 



ci = expOi^) , c 2 = y-y- exp(m 2 t) , (31) 

where ~K\ = iri(t) and 7r 2 = ^(t) are the phases and 

s = \c 2 \ 2 -\c x \ 2 (32) 
is the population difference. We also introduce the phase difference 

X = 7T 2 — TV i + 7 + Aut , (33) 

employ the notation 

a = i (a i2 + a 2 i) , p 12 =pe 11 , (3 = \(3 12 \ , 5 = Au + i (a 12 - a 2 i) , 

and define the dimensionless parameters 

j3 5 



a a 



(34) 
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Then Eqs. (30) can be reorganized to the system of two equations 





= cosx + 5 



(35) 



for the real functions (32) and (33), with the initial conditions s = s(0) and x = x(0). 
Here we have considered the resonant field of the form (22), with the amplitudes depending 
on the spatial variable r but not depending on time t. In general, we could keep in mind 
that the amplitudes Vi(r,t) are slow functions of time, that is, varying much slower that the 
oscillation period 2ti/uj. If so, the value (3 in Eqs. (35) would be a function of time, which 
would permit us to study the behaviour of s(t) and x(t) under the influence of resonant 
pulses. This would enrich the variety of admissible solutions yielding the appearance of 
different transient phenomena, analogously to the case of transient coherent phenomena in 
systems of resonant atoms [53] . 

4 Dynamic Resonant Effects 

The resonant Bose condensate, though being analogous to a resonant atom, possesses, be- 
cause of its collective nonlinear nature, several features that make it rather different from 
the latter. Some of the interesting effects exhibited by the resonant Bose condensate, are 
surveyed below. 

4.1 Mode locking 

We shall say that the modes are locked if the fractional populations |cj(t)| 2 oscillate in time 
so that they do not cross the line |ci| 2 = | C2 1 2 = 1/2. Depending on initial conditions, it may 
be that |ci| 2 > 1/2 and |c 2 | 2 < 1/2, or that |c 2 | 2 > 1/2 while |ci| 2 < 1/2. With regard to the 
population difference (32), the mode locking means that, depending on initial conditions, 
either 



Such a mode- locked regime exists for some region of the parameters (34), in particular, when 
both b and e are small [4]. 

4.2 Critical dynamics 

When the amplitude of the resonant field b or detuning e increase in the magnitude, they may 
reach the values at which the dynamics of the fractional populations experience dramatic 
changes [5,8,12]. This happens on the critical line of the parametric manifold. More precisely, 
for each given set of initial conditions, s , x , there is a separate critical line on the plane 
{b,e}, described by the equation 



1 < s(t) < 



(36) 



or 



< s(t) < 1 . 



(37) 




(38) 
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The critical behaviour of fractional populations, when crossing the critical line (38), has been 
investigated numerically [5,8,12]. The dynamic critical phenomena appear when the motion 
changes from the mode-locked type to mode-unlocked type. This happens when a given 
starting point of a trajectory is crossed by a saddle separatrix on the plane of the variables 
s(t), x(t). Thus, for the starting point s = —1, x = 0, and the detuning e = —0.1, the 
critical pumping amplitude, defined by Eq. (38), is b c = 0.39821. The motion, starting at 
So = —1, %o — 0, is mode locked for b < b c , and becomes mode unlocked for b > b c . The 
phase portrait for b = 0.51 > b c is shown in Fig. 4, where it is seen that the motion beginning 
at «o — — 1, xo — is mode unlocked, i.e. the trajectory lies in the interval — 1 < s < 1. 

Another way of explaining the change from the mode-locked regime to the mode-unlocked 
one can be as follows. Let us introduce the function 

h = 2c\c 2 e i ^ U}t+ ^ = e lx . (39) 

And let us measure time in units of a . Then Eqs. (30) can be presented in the form 

^ = i(s + e)h + ibs , ^ = -b(h- h*) . (40) 
at dt 2 

From the first of these equations, we have 

h = lbs + i — . 

s + e \ dt J 

Substituting this in the second of Eqs. (40) gives 

The latter equation is easily integrated yielding 

(s + e f - b (h + h*) = k , 

with the integration constant 

k = (s + e) 2 - b (ho + h* ) = (s + e) 2 - 2b y/l-s% cosx ■ 
Differentiating the second of Eqs. (40), we get 

^ = -l(s + e)(h + h*)-bs 2 , 

which can be transformed to 

From here, for the shifted population difference 

s = s + e , 



11 



we find 

(§)' = - I^T 1 "* 2 ^- (41) 

with the integration constants 

1 k\ 

fci = fc - 26 2 , A; 2 = Sq + - Sq - y Sq - 2& 2 es , 

where 



so = 2 & (^o - K) = ~b V 1 _ s o sin • 

The exact solution to Eq. (41) can be expressed [54] as a ratio of the Jakobi elliptic functions 
[55]. 

To simplify the consideration, let us set the zero detuning 5 = 0. Then s = s. For the 
integration constants, we get 

ki = si — 2b \jl — Sp cosrr — 2b 2 , k 2 = — s\ — ^ s\ + b 2 (l — Sq) sin 2 x . 
Equation (41) reduces to 

The polynomial in the right-hand side of Eq. (42) has the roots 



s 2 = k x ± ^A; 2 + 4A; 2 , (43) 

which prescribe the admissible region of varying s. For instance, in the case of the initial 
conditions s = ±1, we obtain the roots 

si, 2 = ±1 , s 3 ,4 = ±v / I^4P . (44) 

The left-hand side of Eq. (42) is always nonnegative, hence it should be that 

1 - Ab 2 < s 2 < 1 . (45) 

The analysis of the inequalities (45) tells as that for b 2 < 1/4 the mode- locking regime is 
realized, since the motion is locked in one of the regions 

-1 < s < -(1 -Ab 2 ) ) / . 

b 2 <-) (46) 



1 - 46 2 < s < 1 

depending on initial conditions so = ±1. But for higher b 2 , the motion becomes unlocked, 
with the trajectory wandering in the whole admissible region 

-1 < s < 1 (b 2 > 1) . (47) 

The change of the mode-locked regime to the mode-unlocked one happens at the critical 
value b 2 c = 1/4. 
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4.3 Interference effects 

Since the spatial distribution of two different nonlinear modes y?i(r) and ^(r) is different, 
there appears the interference pattern described by 

p int (r, t) = p(r, t) - pi(r, t) - p 2 (r, t) , (48) 

where 

p(r,t) = \ Cl (t) Vl (v)e' lE ^ h + c 2 (t)^)e- lE2t/h 

p,(r,t) = \c J (t)^(r)\ 2 (j = 1,2) . 

The properties of the pattern (48) were studied in Ref. [12]. 

Different spatial shapes of the nonlinear modes and their different related energies result 
in the appearance of the interference current 

j int (r,t)=i(r,t) - ji(r,t) - j 2 (r,t) , (49) 

in which j(r, t) is the total current in the system, while jj(r, t), with i — 1,2 are the partial 
currents of the corresponding modes [12]. 

4.4 Atomic squeezing 

The evolution equations (30) for the effective two-mode system may be rewritten in the spin 
representation as the equations for the collective spin operators 

N 

S a = >Sj , S± = S x i iSy , 
i=i 

where a = x,y, z. Thence, we may consider atomic squeezing connected with the spin 
operators, because of which it is called spin squeezing [56-58]. Thus, the squeezing of S z 
with respect to S± is described by the squeezing factor 

2A\S Z ) 

Qz = | < S ± > | ' (50) 

in which A 2 (S Z ) =< S 2 > — < S z > 2 . For the resonant Bose condensate, we find 

Q Z = VT^, (51) 

with s being the population difference (32). This shows that for all s ^ 0, the squeezing 
factor Q z < 1, which implies that S z is squeezed with respect to S±. The latter, in physical 
parlance, means the squeezing of atomic states. 
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4.5 Multiparticle entanglement 



To study the possible entanglement for a multiparticle system, one has to consider the 
reduced density matrices. A p-particle reduced density matrix is given by 

Pp (n ...v p ,v[... r' p , t) = Tr ^(n) . . . ^{v p )p{t)^{v' p ) . . . VVi) , (52) 

where rj are Cartesian coordinates, t is time, ip(r) and V^( r ) are field operators, pit) is a 
statistical operator, and the trace is taken over the Fock space. As a basis for calculating the 
trace, we may accept coherent states [48] or Fock-Hartree states [59,60]. At low temperature, 
when practically all N atoms are Bose-condensed, the trace (52) is concentrated on the 
coherent states normalized to the number of particles N. These coherent states have the 
form of the column 

in which 

Vna (r) = VN <p n (r) e ia (0 < a < 2tt) , (54) 

<p n { r ) being the solution to the eigenproblem (1). The phase a here is random, because 
of which the vector (53) may be called the random-phase coherent state or mixed coherent 
state [48,61]. These states are asymptotically orthogonal, 

< VmalVn^ > ^ 5 mn 5 a/3 {N > l) . 

The closed linear envelope 

T = Z{\r] na >} (55) 

forms a subspace of the Fock space. On this subspace, one has the asymptotic, as iV — > oo, 
resolution of unity 

/ IVna >< Vna\ T~ ^ 1 , (56) 
n J 27T 

understood in the weak sense. 

Let us also define a Hilbert space Ti as a closed linear envelope Ti = £{(p n (r)} sup- 
plemented with a scalar product. And let us introduce a p-fold tensor product HP = 
U®U® . . .®H. Then the functions 

p 

V?n(ri,...,I>) EE n^™( r i) 

i=l 

are approximately orthogonal on H p in the following sense. The scalar product 

p 

(VC^D = H(<Pm,<Pn) 
i=l 

has the property 

o< 1(^4,^)1 <i, 



-N/2 



A;=l 



(53) 
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being equal to unity if and only if m — n. But for m ^ n, one has 



Km 1(^,^)1 = (m^n). 

Assuming that the trace in Eq. (52) is concentrated on the condensate subspace (55), 
we find 

Pp (n . . . r p , r[ . . . r' p , t) = £ D n (t) f[ y? n (r,) <P*M , (57) 

n i=l 

where the coefficients are 

D n (t) = < r] na \p(t)\ r] na > . (58) 

The latter can be defined from the normalization conditions 

r TV 1 
J p p (ri . . . r p , ri . . . r p , t) dr ± . . . dr p = — — - , =1, 

which suggests 

D n(t) = j^yKit)? ■ (59) 
For the resonant Bose condensate, with N ^> 1, one gets 

~ iv* , L» 2 (t) ~ i±£ . (60) 

The density matrix (57) cannot be presented as a product of single-particle density ma- 
trices, which implies that the state of the system is entangled [62,63]. This happens for all 
\s\ 7^ 1. Maximal entanglement occurs at s = 0, when Di(t) = D 2 (t) ~ N p /2. 

5 Conclusion 

The method of resonant formation of nonlinear coherent modes makes it possible to create a 
novel type of systems, the resonant Bose condensate. The latter shares many analogies with a 
resonant atom. But, being a collective system, the resonant Bose condensate also displays the 
features essentially distinguishing it from a single resonant atom. First of all, the spectrum 
of the nonlinear coherent modes is described by the stationary Gross-Pitaevskii equation, 
which is a nonlinear Schrodinger equation. These modes, though formally looking as single- 
particle functions, correspond to collective states of a coherent multiatomic system. It would 
be admissible to say that such coherent modes are single- quasiparticle wave functions, but 
not single-particle ones. Here, one should imply a quasiparticle in the sense of Landay, that 
is, an effective dressed single object inside a multiparticle ensemble. 

The evolution equations for the resonant Bose condensate can be reduced to the system 
of two equations for the complex amplitudes of two fractional populations. This reduction to 
an effective two-mode problem is analogous to that happening for a resonant atom. As in the 
case of any two-mode, two-level, or two-component system, the problem allows for the usage 
of quasispin representation. Employing this representation, one could talk about quasispin 
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dynamics, qiasispin waves [64], having some formal properties similar to the dynamics of 
real spins [65], and so on. To study phase portraits, it is also convenient to work in the 
representation employing two real quantities: population difference and phase difference. 

The collective nonlinear nature of the resonant Bose condensate yields the existence of 
several novel effects making this system rather different from a single resonant atom. The 
most interesting of these effects, we have investigated, are: mode locking, critical dynamics, 
interference effects, atomic squeezing, and multiparticle entanglement. These effects can find 
a broad range of various applications. 
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Figure captions 



Fig. 1. Local residual (12) as a function of the dimensionless radial variable r for the 
self-similar approximant (solid line), optimized Gaussian approximation (dotted line), and 
Thomas-Fermi approximation (dashed line) for different values of the coupling parameter: 
(a) g=25; (b) g=250. The residual for the Thomas-Fermi approximation diverges at the 
edge of the atomic cloud. 

Fig. 2. Atomic density (13) as a function of r for the self-similar approximant (solid 
line), optimized Gaussian approximation (dotted line), and Thomas- Fermi approximation 
(dashed line), at the coupling parameter g = 2500. 

Fig. 3. Percentage errors of the self-similar root approximants e\(g) (solid line), e^g) 
(dashed line), and e\{g) (dotted line), as functions of the coupling parameter g, for vortex 
nonlinear modes with n = j = and different winding numbers m and for different trap 
aspect ratio: (a) v — 1, m — 1; (b) v — 1, m — 10; (c) v — 10, m — 0; (d) v — 10, m — 1; 
(e) v = 10,m = 10; (f) v = 100, m = 1; (g) v = 100, m = 10. 

Fig. 4. Phase portrait on the plane of the variables s(t) and x(t) for the detuning e = —0.1 
and the amplitude b = 0.51 that is larger than the critical value b c = 0.39821 associated 
with the initial conditions s — — 1, x — 0. 
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